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Abstract. In shear flows like pipe flow and plane Couette flow there is an extended 
range of parameters where linearly stable laminar flow coexists with a transient 
turbulent dynamics. When increasing the amplitude of a perturbation on top of the 
laminar flow, one notes a a qualitative change in its lifetime, from smoothly varying 
and short one on the laminar side to sensitively dependent on initial conditions and 
long on the turbulent side. The point of transition defines a point on the edge of chaos. 
Since it is defined via the lifetimes, the edge of chaos can also be used in situations when 
the turbulence is not persistent. It then generalises the concept of basin boundaries, 
which separate two coexisting attr actors, to cases where the dynamics on one side 
shows transient chaos and almost all trajectories eventually end up on the other side. 
In this paper we analyse a two-dimensional map which captures many of the features 
identified in laboratory experiments and direct numerical simulations of hydrodynamic 
flows. The analysis of the map shows that different dynamical situations in the edge 
of chaos can be combined with different dynamical situations in the turbulent region. 
Consequently, the model can be used to develop and test further characterisations that 
are also applicable to realistic flows. 
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1. Introduction 

The transition to turbulence in systems like plane Couette flow or pipe flow differs 
from the better understood examples of Taylor- Couette or Rayleigh-Benard flow in 
that turbulent dynamics is observed while the laminar flow is still linearly stable 
(Grossmann 2000, KersweU 2005, Eckhardt et al. 2007, Eckhardt 2008). Evidently, the 
two types of dynamics coexist for the same parameter values. This suggests a subcritical 
transition scenario, where the turbulent state forms around the node in a saddle- 
node bifurcation. Indeed, various bifurcations of saddle-node type have been found 
in these systems (Nagata 1990, Nagata 1997, Clever & Busse 1997, Waleffe 2003, Wang 
et al. 2007, Eckhardt et al. 2002, Faisst & Eckhardt 2003, Wedin & KersweU 2004, Pringle 
& KersweU 2007, Eckhardt et al. 2008) but at least in pipe flow they differ from the 
standard phenomenology in that the node state is not stable but has unstable directions 
as well: it is like a saddle-node bifurcation in an unstable subspace. Numerical studies 
of pipe flow (Schneider et al. 2007) and some simplifled models (Skufca et al. 2006) show 
that also the 'saddle state' has peculiar features. In the higher-dimensional space it 
need not be a fixed point, as in the traditional saddle-node bifurcation scenario, but can 
be dynamically more complicated, i.e., periodic or even chaotic. 

In Couette and pipe fiow the turbulent state forming around the node need 
not be an attractor. Indeed, numerical and experimental evidence indicates that at 
least in the transitional regime the turbulent dynamics is not persistent but transient 
(Brosa 1991, Bottin et al. 1997, Bottin et al. 1998, Moehlis et al. 20046, Faisst & 
Eckhardt 2004, Hof et al. 2006, MuUin & Peixinho 2006a, MuUin & Peixinho 20066, 
Peixinho k MuUin 2006, Peixinho k MuUin 2007, Schneider k Eckhardt 2008a). 
Nevertheless, it is still possible to define a boundary between trajectories directly 
decaying into the laminar state and those first visiting the neighbourhood of the chaotic 
saddle. Trajectories on the turbulent side show a sensitive dependence on initial 
conditions and give rise to rapidly varying lifetimes. This suggested the name "edge 
of chaos" for this boundary (Skufca et al. 2006). In the case of the standard subcritical 
transition scenario, this edge of chaos is given by the saddle state and its stable manifold 
(Ott 2002). There is some evidence that for such a behaviour in plane Couette flow 
(Wang et al. 2007, Schneider et al. 2008). In the case of pipe flow numerical evidence 
suggests that the saddle state is not a single flxed point or a travelling wave, but that 
it rather carries a chaotic dynamics (Schneider et al. 2007). 

In order to explore some of the possibilities in a computationally efficient and 
dynamically transparent manner, we turn to a specifically designed model system. In the 
following we will describe a two-dimensional map that shows much of the phenomenology 
observed in transitional pipe fiow, and at the same time has parameters that allow us 
to discuss the transitions and crossover between different kinds of dynamical behaviour. 
We use the model to study the boundary between laminar and turbulent dynamics, 
and the dynamics in this boundary. In particular, we will argue that the edge of chaos 
and the edge states introduced in Skufca et al. (2006) and Schneider et al. (2007) are 
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the natural extension of the basin boundary concept to situations where the turbulent 
dynamics is transient. 

Studying boundaries of basins of attraction has a long history in dynamical 
systems. It goes back to Cayley for the case of Newton iteration, and to Julia and 
Fatou for dynamical systems defined in the plane of complex numbers (Peitgen & 
Richter 2000, Devaney 2003). To make contact with differential equations much follow 
up work focussed on the conceptually simplest systems of flows in three dimensions, or 
equivalent ly 2d invert ible maps. In principle the generic properties of the boundaries 
between the domains of attraction of different types of invariant sets (sinks, saddles 
nodes, limit cycles, and chaotic sets) have exhaustively been classifled (Robert et al. 
2000, Ott 2002) for these systems by considering (i) the possible sections of the respective 
stable and unstable manifolds and (ii) the possible impact of (dis-)appearance of stable 
orbits in saddle-node bifurcations. However, careful inspections of the parameter 
dependence of 'explosions', where the features of invariant sets alter qualitatively, 
can occasionally still unearth surprises in systems as simple as the Henon map 
(Osinga 2006). Higher-dimensional chaos ( "hyperchaos" ) shares common themes with 
low- dimensional chaos (Rossler 1983), but there also are important differences due to 
the additional freedom of changing dynamical connections between chaotic sets (Grebogi 
et al. 19836, Lai & Winslow 1995, Dellnitz et al. 1995, Ashwin et al. 1996, Kapitaniak 
et al. 2003, Rempel et al. 2004, Pazo & Matias 2005, Tel & Lai 2008). Besides fluid 
mechanics other important flelds of applications of hyperchaos are transition state theory 
(Kovacs & Wiesenfeld 2001, Wiggins et al. 2001, Waalkens et al. 2004, Benet et al. 2005) 
and the quest for the (domain of) stability of irregular and stable synchronised states in 
systems of coupled oscillators (see Pikovsky et al. (2001) for an overview). Considerable 
insight in the latter problem come from studies of two symmetrically coupled logistic 
maps (Yamada & Fujisaka 1983, Fujisaka & Yamada 1983, Gu et al. 1984, Pikovsky & 
Grassberger 1991, Maistrenko et al. 1998, Kapitaniak et al. 1999, Kapitaniak et al. 2003). 
More recently also the generalisations to asymmetric coupling (Hu & Yang 2002, Kim 
et al. 2003) and more complex maps (Lai 2001, Kim et al. 2003, Ashwin et al. 2004) 
have been explored. 

The present study is motivated by observations on the turbulence transition in 
situations where the laminar proflle is linearly stable, and hence will use descriptions 
like 'laminar' and 'turbulent' to describe the two dominant state between which we would 
like to determine the basin boundary or edge of chaos. One of our principle interests 
will be in situations where the dynamics on the edge of chaos separating (transient) 
turbulence and laminar motion is chaotic. To that end our model must have at least 
two continuous degrees of freedom — one degree of freedom for the dynamics in the 
edge, and a second one perpendicular to it. A minimal model of the phase-space flow 
would then require at least a four-dimensional invert ible map, but then we would loose 
the advantages of the graphical representation of the invariant sets and their domains 
of attraction that are available in lower dimensions. As in the approaches to model 
synchronisation of coupled nonlinear oscillators, we will therefore design a system of 
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two coupled 1-d maps. 

The paper has three main parts. In the first part (section[2]) we introduce the model, 
discuss the dynamics of the uncoupled case, and introduce the considered coupling. The 
second part (sections [3] and ^ deals with the dynamics of two coexisting attractors: 
In section |3] we discuss the shape and dynamics of the attractors, and the transient 
dynamics in the respective basins of attraction. Section [4] addresses the dynamics of the 
relative attractor on the basin boundary between the attractors, and how this dynamics 
effects the shape of the separating boundary. In the third part of the paper we turn 
to the case of a chaotic repellor in the turbulent dynamics which mimics turbulent 
transients decaying to a laminar flow proflle: Section [5] deals with the case of a chaotic 
saddle coexisting with a flxed point attractor. We discuss the metamorphosis of the 
basin boundary at the crisis where the attractor turns into a chaotic saddle. Finally, in 
section [6] we conclude the paper with summarising remarks and discuss how the flndings 
on this 2D model relate to observations in shear flows such as turbulent pipe and plane 
Couette flow. 



2. The two-dimensional map 

To admit coexistence of laminar and a turbulent dynamics one degree of freedom of 
the map must be chosen along a phase-space direction separating regions with these 
different types of dynamics. A second degree of freedom is needed to capture the 
dynamics perpendicular to this direction, and to allow for dynamics within the boundary 
between laminar and turbulent dynamics. We think of the two coordinates of the map 
as representing the energy content of the perturbation (x-direction) and the dynamics in 
an energy shell (y-coordinate). The x-coordinate interpolates between a laminar and a 
turbulent dynamics. The ^-coordinate models all other degrees of freedom. In the latter 
direction the map is globally attracting towards a region near y = I. The combined map 
has a flxed point, corresponding to the laminar proflle, and — for suitable parameter 
values — also a region with a chaotic dynamics corresponding to turbulent behaviour. 
In the following we flrst describe the two uncoupled maps in x and ^, and then we 
discuss their coupling and its consequences for the dynamics. 



2.1. Dynamics in x 

For the dynamics along the energy axis, we use a map that has a stable flxed point at 
X = —2, and a chaotic dynamics for x > 0. The former corresponds to laminar flow, 
and the latter mimics turbulent motion. An intermediate flxed point at x = separates 
the laminar region x < from the turbulent one at x > 0. It is unstable. These features 
are contained in the one-parameter map [flgure[l](a)] 

Xn+i = f{xn]a) (la) 
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Figure 1. The map along the x-direction. (a) The red hne shows the function /(x; a) 
for a = 3.8, and the green ones indicate the evolution of two trajectories starting at 
X = ±0.1, respectively. The maximum of the map is always at x = 1/2, irrespective 
of the value of a, and takes the value f{x = 1/2; a) = a/4. The initial condition 
0.1 approaches the chaotic attractor in [/(1/2); /^(l/2)] = [0.18; 0.95], and the one 
starting at —0.1 approaches the fixed point at x = —2. (b) Bifurcation diagram of the 
map f{x;a). Green dots represent points on either of the two attractors of the map, 
and the red lines unstable fixed points. 



with 

, \ axil - x) for X > X* = ( 1 - a/1 + 8/a] / 2 , 

f{x]a) = < V ^ // (16) 

[ -2 else 

Here x^ is the leftmost intersection between the constant value —2 for x < x*, and the 
quadratic part at x > x*. With this choice the map is continuous. 

The bifurcation diagram for this map is shown in figure [l](b). We will only be 
interested in parameter values a > 1/3 where > — 2. In this case the map has a 
stable fixed point at x = — 2, which absorbs all initial conditions starting outside the 
interval [0, 1]. Over the interval x G [0, 1] the map coincides with the logistic map and 
shows its familiar bifurcation diagram. For all 1/3 < a < 1 there are stable fixed points 
at X = — 2 and x = 0. In addition, there is an unstable fixed point at x^ = 1 — 1/a, 
which lies between —2 and 0. At a = 1 the fixed point Xg crosses x = 0, and the 
two fixed points change stability in a transcritical bifurcation. For a > 1 the point 
X = is unstable, and Xg is a stable fixed point. At a = 3 the fixed point Xg undergoes 
a first period doubling, and subsequently follows the period-doubling route to chaos. 
Beyond a 3.59 there are chaotic bands extending from /(l/2;a) = a/4 down towards 
f{l/2;a) = f{a/4;a) = {a/2f{4-a). 

At a = 4 the chaotic band generated by the period doubling collides with the 
unstable fixed point at x = 0, leading to a boundary crisis (Grebogi et al. 1982, Grebogi 
et al. 1983a, Grebogi et al. 1987, Ott 2002). For a > 4 some points near the maximum 
of the parabola are mapped outside the interval [0, 1] and the attractor turns into a 
chaotic saddle. All points except for a Cantor set of measure zero will eventually map 
outside the interval and then be attracted to the laminar fixed point at x = — 2. The 
Cantor set contains an infinity of orbits which follow a chaotic dynamics and never leave 
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Figure 2. The map in the ^-direction, (a) The function g{y] b) for 6 = 6, and 
the trajectory of y = 0. For ah b the maximum of g{y;b) is at {y^g{l;b)) = (1;2). 
Moreover, g{2; b) = g{0; b) = 2/(1 + b) such that the attractor of the map always hes 
in the interval [2/(1 + 6); 2]. (b) Bifurcation diagram of the map g{y] b). As in Fig. [l] 
green dots represent points on the attractor, and the red lines mark two unstable fixed 
points of particular interest. 



the interval (c/. Tel 1990, Ott 2002). 

In summary, depending on the parameter values, the x-map shows the coexistence 
of a stable laminar state with one of three possible types of non-laminar dynamics: 
another fixed point, a chaotic attractor, or a chaotic saddle. The coexistence of a stable 
laminar fixed point at x = — 2 with a transient chaotic dynamics in the map for a > 4 
mimics the coexistence of a transient turbulent dynamics with a linearly stable laminar 
steady flow. The direct domain of attraction of the laminar state at x = — 2 is bounded 
towards positive x by an unstable flxed point at x = 0. 



2.2. Dynamics in y 

The ^/-dynamics represents the motion within the energy shell. In the simplest case it 
is globally attracting towards a globally stable flxed point. Then only the x-dynamics 
matters, and it represents the dynamics along its unstable direction. In order to model 
the motion in the energy shell we consider a unimodal {i.e., a single-humped) map of 
Lorentzian type (flgure[2](a)) that maps large \y\ towards the region ?/ — 1, 

Vn+i = giVn'M , (2a) 

with 

In its flrst iteration the map collects all initial conditions into the interval [0,2]. In 
this interval the map can have up to three flxed points y^. For the discussion of the 
properties of the map and the flxed points, it is convenient to solve the flxed point 
equation for the parameter and to study 
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By evaluating dbo/di/p = one verifies that there is a saddle-node bifurcation at the 
critical value i/cr = (3 — a/5) /2 0.382. This corresponds to the parameter value 
bcr = bod/cr) — 11.09. Consequently, there is only a single fixed point for b < bcr-, and 
there are three fixed points for larger values of b. 

Making use of equation ^ in order to evaluate dg{y]b)/dy = —1 one verifies 
that the single fixed point is stable for y > ys = (l + a/5) /2 1.6182, i.e., for 
b < bg = bo{ys) — 0.618. Beyond bg the fixed point undergoes a period-doubling route 
into chaos, and produces a broad chaotic band in the interval [0,2]. At bcr there is a 
saddle-node bifurcation in the support of the attractor, which transforms the attractor 
into a saddle. For larger values of b this saddle coexists with a globally stable fixed 
point. 

For later reference we introduce also the Lyapunov number A of the map, which 
describes how a small distance 6yo = |^o^^ — ^o^^| between two close-by initial conditions 
^0^^ and ^/o^^ grows with the number j of iterations, 

5y^ = \yf-yf\^SyoA'. (4) 

The Lyapunov numbers can be defined for invariant sets, such as the maximal chaotic 
invariant set (Ac) and for the attractor (A^). The distinction is important whenever the 
two numbers do not coincide, as in cases where an attracting periodic orbit is surrounded 
by an invariant chaotic set. The two Lyapunov numbers for the map (26) are shown in 
figure [3l The Lyapunov number for the maximal chaotic invariant set is shown as a solid 
red line: it always remains above 1. The Lyapunov number of the attractor is shown by 
a dotted green line. It takes values smaller than unity in the parameter windows where 
there is an attracting periodic orbit. 

In summary, the main features of the ^-dynamics are that it is globally contracting 
towards the interval [0;2], and that depending on the parameter values one can have 
one of three types of invariant sets: (i) a stable periodic orbit of period 2^ with n = 
(i.e., a fixed point) for b < bg and larger n in the subsequent period doubling cascade; 
(ii) a chaotic attractor for numerous parameters in the range bs < b < bcr] or (iii) a 
chaotic saddle coexisting with a periodic orbit (in the periodic windows of the previous 
parameter regime) or a fixed point for bcr < b. 



2.3. The coupling 

Without a coupling between the two maps, the three possibilities in the x-dynamics 
combine with the three possibilities in the ^-dynamics for nine different regimes. Now 
we introduce a coupling between both degrees of freedom. The specific form of the 
coupling should not be important if it preserves a few properties. For instance, we want 
to keep a locally stable fixed point for the laminar state also in the coupled dynamics. 
Specifically, the y-map should have a stable fixed point at x = —2. We therefore 
introduce an x-dependence in the parameter b of the ^-map such that the coupling 
vanishes for x — 2, thereby maintaining the stability properties of the uncoupled map: 

6(x) =7(2 + x). (5) 
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(a) 




(b) 




Figure 3. Lyapunov number A of the attracting set (green) and the chaotic invariant 



set (red) of the map (2h). The former number is obtained from the last 2-10^ iterates of 
a trajectories that is 10^ iterations long, and the latter by means of the thermodynamic 
formalism [cf. Tel (1988, 1990) for a description of the algorithm], (b) Magnification of 
a small parameter interval to show that in regions where, within numerical uncertainty, 
there are no stable periodic orbits, the two Lyapunov numbers coincide. 



We refer to this fixed point as the laminar fixed point. Since the non-trivial x-dynamics 
lies within the interval [0, 1], the range of h values varies between 27 and 87, so that the 
parameter 7 selects the type of ^—dynamics for the chaotic regime in the x dynamics. 

To complete the coupling we also introduce an influence of the ^-dynamics on the 
x-dynamics, since otherwise the bifurcations are determined by the x-map alone: we 
shift Xn by the deviation of yn from the position of the maximum before applying the 
mapping, i.e., 

Xn+i = f{xn - e{yn-l)] a) (6a) 
Vn+i = g{yn]b{^n)) (66) 



with the speciflc forms (16), (26), and (??) for /(x), g{x)^ and 6(x), respectively. In 
this paper we will concentrate on the weak-coupling limit where e <C 1. Unless stated 
otherwise this parameter will always take the value e = 0.03. 

This completes our deflnition of the coupled map. Through appropriate choices 
of the parameters a and 7 we can — one by one — study the nine parameter regimes 
with their qualitatively different dynamics. We here begin with the six cases where 
the non-laminar x-dynamics is attracting, and a laminar and a non-laminar attractor 
coexist. The case of a transient dynamics will be taken up in section [5| 

3. Two coexisting attractors 

3.1. Attractors and basins 

Figure |4] shows the domain of attraction of the laminar flxed point together with the 
non-laminar attractor. The panels on the left-hand side refer to a = 1.2 immediately 
beyond the crossing of stability, where /(x,a) has a stable flxed-point attractor at 
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Figure 4. For initial conditions (x, y) the colours indicate the number of iterations 
required to reach the laminar fixed point. The colour code runs from one iteration 
(blue) to more than ten iterations (red). Initial conditions in the white region are 
attracted to the chaotic attract or, which is also shown by red dots. The panels in 
different rows refer to different values of 7: (a,d) 7 = 0.2, (b,e) 7 = 3, and (c,f) 7 = 6, 
respectively. In these cases the ^/-dynamics shows a single fixed point, chaos, and a 
fixed point coexisting with a chaotic saddle, respectively. The left panels (a-c) and 
right panels (d-f) refer to a = 1.2 and a = 3.8, respectively. For the panels (a-c) the 
non-laminar x-dynamics amounts to a fixed point, and for (d-f) it is chaotic. In all 
panels e is set to 0.03, except for the top left one, where e = 0.01, because in the latter 
case the non-trivial fixed point at x > disappears for e > 0.012. 
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Figure 5. Magnifications of the boundary between the region of initial conditions 
approaching the laminar fixed point at (—2,2), and the chaotic attractor at x > 
for 7 = 3, e = 0.03, and different values of a: (a) a = 1.2, (b) a = 1.6, (c) a = 2.0, 
(d) a = 3.8. All initial conditions which approach the attractor to the right are shaded 
in red, and the colour gradient from blue to yellow indicates the number of steps it 
takes to reach the laminar fixed point with the colour coding also used in figure [4j The 
cross-over from a rough to a smooth boundary arises at a :^ 1.6. 



X = 1 — 1/1.2 :^ 0.17. The panels on the right-hand side refer to a = 3.8. In this case 
the fixed point at x = — 2 coexists with a chaotic attractor. Without coupling, for (6a) 
with 6 = the attractor is located in the interval [0.18; 0.95]. 

For 7 = 0.2 [figure |4](a,d)] the iteration of y directly approaches the fixed point 
dit y 1.8, and subsequently only wiggles around this point due to the perturbation 
arising from the x-dynamics. 

For 7 = 6 [figure [4](c,f)] the iteration of y approaches the fixed point at y 0.38. 
However, in this case the fixed point is surrounded by a chaotic saddle, and the approach 
may involve long chaotic transients. 

Finally, for 7 = 3 [figures |4](b,e)] the parameter b{x) varies between 6 and 9 for x- 
values in the interval [0, 1]. For parameter values a slightly below 4 one hence obtains a 
chaotic dynamics for both x and y. In this case the attractor varies over a considerable 
range of ^-coordinates, and one can clearly see the strong influence of the coupling. 
When there are broad chaotic bands in both directions [figure [4](e)] the attractor can 
even extend to negative values of x. 
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3.2. The boundary between the two attractors 

We now focus on the boundary separating the basins of attraction of the laminar and 
the chaotic attract or to the left and right, respectively. If e = 7 = 0, it coincides with 
the y-axis: All initial conditions with x > are attracted to the turbulent dynamics, 
and the ones with x < to the laminar state. Moreover, all points with x = are 
immediately mapped into the hyperbolic fixed point (0,2). The hyperbolic fixed point 
then becomes a relative attractor, since it is an attractor for initial conditions in the 
boundary between the two attractors. For e and 7 nonzero but small, the hyperbolic 
point is slightly shifted, and the boundary no longer coincides with the y-axis, but it 
remains smooth. The boundary can be determined by picking initial conditions with, 
say, prescribed y and varying x and following them for some iterations forward in time: 
it can then be bracketed by a pair of x- values where one initial condition iterates towards 
the laminar state and the other towards the turbulent one. This method allows us to 
track the dynamics in the boundary not only in the case where the relative attractor 
is a fixed point, but also when it is more complicated. In a hydrodynamic setting this 
approach has been explored in the framework of low-dimensional shear fiow models 
(Skufca et al. 2006) and direct numerical simulation of pipe fiow (Schneider et al. 2007). 

In figure |4] the boundary between the two attractors is the boundary of the region 
shaded from blue to yellow. It appears to be smooth for a = 3.8 and for 7 = 0.2. In 
contrast, it looks irregular for a = 1.2 and 7 = 3 or 6 [figures [4](b,c)]. The magnifications 
in figure [5] confirm the roughness of the boundary and indicate a crossover from a smooth 
to an irregular boundary as a decreases from 3.8 to 1.2, with 7 = 3 and e = 0.03 fixed. 

There are two elements needed to understand the emerging roughness of the 
boundary: the first one is the observation that states in the boundary are attracted 
to a subset of the boundary itself, i.e., the dynamics in the basin boundary converges 
to an edge state. The second observation is that when the edge state is chaotic a rough 
boundary can form provided that the Lyapunov exponent for the chaotic motion on 
the basin boundary is larger than the one characterising the escape from the boundary. 
These two aspects are discussed in the next section. 



4. Edge states and relative attractors 
^.7. Identifying the edge state 

In order to follow a trajectory for long times and to be able to identify the relative 



attractor, the bracketing of trajectories in the boundary described in section |3.2| has 
to be refined after some time. After all, the distance between the trajectories in the 
pair bracketing the trajectory on the boundary grows exponentially with the number 
of iterations. Specifically, we proceed as follows. We take initial conditions for the two 
trajectories that have equal ^-values and x- values separated by less than 10~^. The 
two trajectories are followed until max(5xj,5^j) exceeds 5 • 10~^. Then a new pair is 
determined with yo = [y^ + y^)/2 and Sx^ < 10~^. An alternative approach could 
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start from the observation that the hne connecting the two trajectories wiU be oriented 
along the direction of the largest Lyapunov exponent of the map and search for a 
refinement along this line. Here and in the previous applications to pipe flow (Schneider 
et al. 2007, Schneider & Eckhardt 20086) it was observed that the flrst approach, which 
repeatedly projects the line segment between the two points to a flxed direction in space, 
is more robust and converges more reliably, especially in cases where the geometry of 
the boundary is complex. 

The dynamics in the edge state is explored further in flgure [6| Presented are 
two situations where the boundary shown in flgure [5] appears smooth (a = 1.2, left 
column) and rough (a = 3.8, right column), respectively. The two frames in flgure |6](a) 
show trajectories on the boundary constructed by the edge tracking algorithm. The 
trajectories nicely reproduce the features of the boundaries also shown in flgure [5](a) 
and (d). The difference between the two flgures is that the boundary in flgure [5] emerges 
from a two-dimensional search, whereas the one in|6]is determined by a following a single 
trajectory. This allows us to show the time series of the coordinates of edge trajectories 
and the associated return maps for the ^-coordinates in row (b) and (c), respectively. 
By visual inspection it is very hard to see differences to the unperturbed dynamics of 

To demonstrate effects introduced by the coupling of the dynamics to the unstable 
x-direction we subtract the functional form of the ^-map. The deviations from the 
unperturbed y-dynamics, 5yn = Vn+i — g{yn)^ differ substantially for smooth and rough 
boundaries: For a = 3.8 the iterates lie on a smooth, double valued curve. Its double- 
valuedness reflects the influence of a non-trivial dynamics in x, which follows iterates of a 
map with a single bump, see the iterates in row (e). However, the relation between x and 
y is single valued, and therefore there is not much disorder. For a = 1.2, the distribution 
of iterates looks rather noisy (e), and no simple relation between their images can be 
found. 

Note that in both cases the dynamics in y is chaotic, and along the x-direction close- 
by trajectories escape exponentially from the vicinity of the boundary: both Lyapunov 
numbers are positive. On the other hand, for a = 3.8, the different branches of the return 
map come to lie on a smooth invariant set, while for smaller a the basin boundary is 
a rough invariant set. In the next subsection we argue that this difference is due to a 
crossover of the absolute values of the Lyapunov numbers, just as it has been discussed in 
the context of unstable-unstable pair bifurcations (Grebogi et al. 19836, Tel & Lai 2008). 

4^2. Transition between smooth and rough boundaries 

Close to a = 1.6 the boundary crosses over from a highly irregular geometry to a line 
with only few kinks, whose number progressively decreases for even larger values of a. 
Similar transitions between rough and smooth boundaries have previously been seen 
in unsteady-unsteady pair bifurcations (Grebogi et al. 19836, Ott 2002) and phase- 
synchronised chaos (Hunt et al. 1997, Rosa and Ott 1999). They are related to a 



Vollmer et al — Basin boundary^ edge of chaos, and edge state in a 2d model 13 
(a) 





-0.06 -0.03 



-0.06 -0.03 





Figure 6. Evolution of trajectories on the boundary separating convergence to 
laminar and to turbulent motion (7 = 3 and e = 0.03 for all panels; left panels: 
a = 1.2; right panels: a = 3.8). (a) The trajectories on the boundary trace out the 
shape of the boundary, which is rough for a = 1.2 and smooth for a = 3.8. For later 
reference the right panel also shows the edge of chaos for a = 4.2 beyond the crisis of 
the attractor. It has been shifted by {x;y) = (0.001; —0.04) to the lower right, (b) 
When started at x = the trajectories rapidly converge to a chaotic attractor located 
in the interval y G [2/7; 2]. (c) Return map for the ^/-coordinate. Within numerical 
accuracy it agrees with g{y;6) (solid red line). The first ten iterations are explicitly 
indicated — the points visited during the initial 10000 time steps are indicated by green 
crosses (+). (d) Deviation of the values plotted in (c) from the function g{y;6). (e) 
Return map for the x coordinate. The first iterations and subsequent points indicated 
as in panel (c). 
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crossover of the two Lyapunov numbers of the map. To gain insight into the transition we 



estimate the slope |5^o/^^o| of the boundary at a point (xq, ?/o)- Linearising equation (6a) 
around the iterates (xj^i/j) of the considered point we find 

Sxj+i = f {xj - e{yj - 1)) {5xj - eSyj) = aj {5xj - eSyj) (7) 

where a = f'{x) = a(l — 2x) is the derivate of f{x]a). Since aU points he on the 
boundary, Xj — e {yj — 1) is close to zero for all as can also be verified by inspection 
of figures [i] and [5| Consequently, is always evaluated at a point close to zero, and aj 
takes values close to a. By recursively working out equation ([T]) we find 

SXn= I (Jn-k 1 Sxq - e ^ I an-k 1 ^Vn-j 

\k=l J j=l \k=l J 



« = + e 2 (8) 

1 h=l ^n-k I lk=l ^n-k 

For an initial perturbation which is located on the boundary the deviation Sx^ is 
bounded, and — in the present case — in absolute value it is much smaller that unity, 
(c/. figure [5]). On the other hand, for large n the denominator H^^i ^n-k takes on very 
large values — after all, a > 1 and cjk — Consequently, 

n-j-l ^ 

S^o^e V (9) 



=0 



In the limit of very small perturbations and large n we can approximate the product 
by its asymptotic scaling, i.e., 

n-l 
k=0 

In addition, according to equation ^ the parameter b of g always takes values very 
close to 27 because all Xj are very close to zero. As shown in figure [6](c) the dynamics 
of the y coordinate essentially amounts to the unperturbed dynamics such that we may 
use equation ^ to related Syj to Syo. 

In the scaling regime the sum in equation ^ can be worked out, yielding 



Sxq 



TTn 1 A 

eH— with H=-. (10) 

H — 1 a 



In the limit n ^ oc the right hand side of equation (10) remains finite only if H < 1. 



Hence, the boundary will be smooth for i7 < 1, or A < a. On the other hand the 



bound (10) diverges for H > 1. In this case the slope diverges at least for some points 



on the boundary, which will hence be rough{|] 



X A discussion of the abundance and distribution of singular points, and the fractal dimension of the 
basin boundary lies beyond the scope of the present manuscript. They can explicitly be worked out 
along the lines indicated in Rosa and Ott (1999). 
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As noted above for points on the boundary the parameter h oi g always takes values 
very close to 27. According to figure |3] one thus finds that A c:^ 1.59 for 7 = 3. The 
crossover from a rough to a smooth boundary should therefore occur at a 1.59, which 
is in excellent agreement with the numerical findings of figure [5| 

This completes the characterisation of the attractors and their basin boundary. In 
the following section we address the case of a chaotic saddle coexisting with a laminar 
fixed point. 

5. Transient chaos 
5.1. Lifetime Plots 

The six cases discussed in the preceding sections cover the cases of coexisting attractors. 
However, close to the transition in plane Couette fiow and pipe fiow the turbulent 
dynamics is transient, so that also the cases of a coexistence between a laminar fixed 
point and a chaotic saddle that supports transient chaotic dynamics are of interest. Our 
map realizes this for e = 0.03 and a > 4 (see figure ig. As in fi gure [i] we consider the 
three cases (a) 7 = 0.2, (b) 7 = 3 and (c) 7 = 6. 

When the parameter a exceeds a critical value acr(7), the laminar fixed point 
becomes globally attracting except for a measure zero set containing periodic and 
aperiodic trapped orbits left over from the attractor. This is apparent in the plots 
in figure [7| which show the lifetime of initial conditions (x,^) for a = 4.0 and different 
values of 7. For 7 = 0.2 [figure[7](a)] the critical value is larger than 4.0, i.e., there still 
is a stable chaotic attractor coexisting with the laminar fixed point. However, we already 
see two 'fingers' approaching the attractor from the top and from the bottom. When 
increasing either 7 or a these fingers are joined by additional narrower fingers which all 
simultaneously collide with the attractor at the parameter value (7). Beyond this 
crisis most of the trajectories of the former attractor escape through the regions where 
the collision took place (Ott 2002). The orbits of the attractor which never enter the 
regions form a chaotic saddle. 

The panels figure [7|(b,c) show the situation beyond the crisis. The blue areas iterate 
to the laminar fixed point in one and two iterates, respectively. The dark green strips 
near x and x 1 arrive at the fixed point in three iterations, and initial conditions 
in the widest fingers (also dark green) pointing towards (x,^) = (0.5; 1) escape to 
the laminar fixed point in four iterations. On the next level there are four lighter 
green fingers lying between the widest fingers and the outer regions (0 < x < 0.5 and 
0.5 < X < 1), respectively, which are mapped to the fingers near x 0.5. With 
each additional iteration, the number of fingers doubles. At the crisis all fingers 
simultaneously collide with points lying at the upper and lower boundaries of the 
attractor. They can be interpreted as a primary collision of the attractor with its 
basin boundary, and the simultaneous collision of all the pre-images of this point. 

What happens to the basin boundary of the attractor when going through the 
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crisis? The chaotic attractor embedded in the basin boundary merges with the attractor. 
We have seen that this generates a fractal set of "holes" (actually the fingers) through 
which trajectories of the former attractor escape to the laminar fixed point. The chaotic 
trajectories that never enter the fingers form a Cantor set. Since trajectories starting in 
the domain of attraction are attracted towards (a small neighbourhood of) the Cantor set 
and those starting in the vicinity of this set escape almost certainly to the laminar state, 
the Cantor set forms a chaotic saddle for the dynamics. There are orbits approaching 
this set from outside, but randomly selected points in the vicinity of every point of the 
Cantor set eventually approach the laminar state with probability one. 

Figure [6](a) shows orbits on the boundary separating the respective domains of 
attraction towards the laminar fixed point and the chaotic set. As demonstrated in 
figure [6](a, right panel) these orbits change smoothly when the system undergoes crisis. 
The transition from a system with a chaotic attractor to one with only chaotic transients 
is solely refiected in the fact that the orbits on the edge of chaos attain new pre-images. 
Their forward dynamics is not affected. In this respect the trajectories forming the 
basin boundary remain a well-defined set also beyond crisis. Their closure is the edge 
of chaos. 

Most initial conditions from the former attractor sooner or later cross the edge of 
chaos. On the other hand the close- by points on the Cantor set, which forms the chaotic 
saddle, never cross the edge of chaos. Some of them step on the edge and are attracted 
towards the relative attractor on the edge of chaos. They give rise to the additional pre- 
images mentioned above. Most points of the Cantor set, however, only closely approach 
the edge of chaos, and subsequently follow its unstable directions to explore the full 
support of the Cantor set. In this sense the edge of chaos remains a well-defined object 
also after the crisis. It separates initial conditions where all orbits immediately decay 
to the laminar fixed point from a region where they can perform a chaotic transient — 
either short but occasionally also very long. In this sense the edge of chaos separates 
initial conditions which are characterised by their diflFerent finite-time dynamics rather 
than by their asymptotic behaviour: the notion of the edge of chaos extends the concept 
of a basin boundaries between two attractors to the situation of an attractor coexisting 
with a chaotic saddle. 

5.2. Parameter dependence of the lifetime for initial conditions on the y axis 

A useful and experimentally accessible indicator for the boundaries and their dynamics 
are lifetimes of perturbations, figure [7| shows the lifetimes for fixed parameters and a 
two-dimensional domain of different initial conditions. The frequently used lifetime plots 
for turbulence transitions differ from this one in that they usually show the deviations 
for a combination of one coordinate (the amplitude of a velocity field) and a parameter 
(the Reynolds number). 

To gain insight into the relation between these two kinds of lifetime plots we first 
consider the conceptually simplest case where the lifetime of trajectories starting on the 
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Figure 7. The change of the structure of the invariant chaotic set (red dots) when 
the attractor undergoes crisis. The colour coding indicates initial conditions arriving 
at the laminar fixed point in at most 10 iterations. The tenth iteration of points which 
have not yet reached the laminar fixed point after 40 iterations are indicated by red 
boxes. Parameter values are a = 4.0 and e = 0.03 for all panels, while 7 takes on 
different values: (a) 7 = 0.2 immediately before the transition from the attractor to 
the saddle, (b) 7 = 3, and (c) 7 = 6. 

y axis is plotted as a function of a and y (figure [s]) . The large blue domain in the upper 
half indicates parameters and initial conditions that are quickly attracted to the fixed 
point. The large red region in the lower left indicates initial conditions which never get 
to the laminar fixed point, since the turbulent domain is an attractor. 

The magnification figure |8](b) focusses on the fuzzy regions in the lifetime plot for 
a 1.25. As we have seen in figure[5]the boundary between the two coexisting attractors 
in the coordinate space (x, y) is rough for these parameters. As a consequence the y- 
axis repeatedly crosses the boundary between the domains of attraction of the respective 
attractors. This gives rise to the observed spiky structure of the interface in the a-y 
plot figure |8](b). Beyond a 1.6 the basin boundary is smooth [figure [5](c,d)], and also 
in an a-^-plot there is a sharp boundary between the two domains. It is located close 
toy^ 0.82. 

When the attractor undergoes the boundary crisis at = 3.93 the fingers from 
figure [7] are visible also in the a-^-plot. They form a hierarchical structure of regions 
that are mapped into the crisis region and subsequently rapidly approach the laminar 
state. Note that, when sufl&ciently resolved, also in this case all fingers extend to the 
critical parameter value — 3.93. 

5.3. Generic parameter- coordinate dependence of the lifetime 

In figure |8] we chose a section aligned almost parallel to the edge of chaos. On the other 
hand, in the applications (Darbyshire & Mullin 1995, Skufca et al. 2006, Schneider 
et al. 2007) the amplitude of a perturbation of the laminar state is varied, i.e., initial 
conditions are chosen along a line extending from the laminar fixed point towards the 
phase-space domain admitting chaotic motion. Such a line intersects the boundary more 
or less perpendicularly. In that case one encounters a sharp, smoothly varying boundary 
between the laminar and turbulent regions for all values below the crisis: It is no longer 
possible to resolve the roughness of the boundary close to a 1.25. In view of this we 
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Figure 8. Metamorphosis of the x = section through the boundary for fixed 7 = 3.0, 
e = 0.03, and varying a. Panel (a) gives an overview, (b) a blow-up of the region with 
a rough basin boundary, and (c) focusses on the transition of the chaotic attractor to a 
chaotic saddle. The colour coding is the same as in figure [5] vertical sections through 
the present plots exactly agree with a section along the y axis of the corresponding 
plot in figure [5] 



focus on the region close to the crisis. The appropriate parts of the parameter plots for 
four different slopes m of the line 

?/ = 2 + m (x + 2) (11) 

are shown in figure [9j 

All panels of figure [9] show hierarchical organised traces of the fingers that we 
also saw in figure [8](c). This shows that folded and hierarchically organised structures 
in lifetime-plots are generic. They do not dependent on the specific direction along 
which initial conditions are chosen. On the other hand the choices differ in the detailed 
structure of the folds: Figure [9](a) shows the situation where the initial conditions 
on the line approach the saddle, but do not intersect it. In this case the folds are 
nicely aligned, and they extend down to different parameter values a well below the 
bifurcation. After all [cf. figure [7|(a)], the fingers invade the domain of attraction before 
they collide with the attractor at the crisis, and the tips of the finer fingers come down 
at a later time. Figure [9](b) corresponds to the situation where the line touches the 
outer edge of the saddle. Consequently, it is exactly along this line that all finger tips 
simultaneously collide with the attractor. Before the crisis, all initial conditions proceed 
into the attractor, and at the crisis there is a fractal set of folds with initial conditions 
escaping to the laminar state appearing all at once. Subsequently, only the scaling of the 
width of the folds, and hence the fractal dimension of the remaining saddle changes. In 
figure |9](c) the initial conditions giving rise to chaotic motion lie right in the heart of the 
chaotic invariant set. In this case the folds also appear simultaneously at the crisis. A 
new feature is that the internal dynamics of the saddle gives rise to a non-trivial bending 
of the folds. For many values of x, in particular x = 0.5, there is not a unique value of 
a which separates regions of persistent chaotic motion for smaller a from a decay to the 
laminar state. Rather, there can be multiple switching between these possibilities as a 



is increased. When the line (11) intersects the chaotic set only at its lower boundary 



[figure [9](d)] the qualitative features of the position-parameter plot are the same as in 
case (c), except that the multiple switching is less pronounced. 
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Figure 9. Metamorphosis of a section through the boundary where the coordinates 
are varied along a hue ( pTj ). The different panels correspond to (a) m = 1/3, (b) m = 0, 
(c) m = — 1/3, and (d) m = —2/3, respectively. All other parameters and the colour 
coding are the same as in figure [8| 

In all cases the observed structure of folded hierarchical tongues are reminiscent of 
the observations in studies of minimal perturbation amplitudes in pipe flow (Darbyshire 
& Mullin 1995, Schneider et al. 2007). Thus, this model provides further support for 
the idea that transient turbulent motion is generated by a chaotic saddle that coexists 
with a laminar fixed point in the state space of linearly stable shear flows. The following 
section discusses in more detail the implications of these findings to the transition to 
turbulence in linearly stable shear flows. 

6. Discussion 
6.1. Methods 

We have suggested a low-dimensional model in which we can analyse methods and 
concepts that recently been used in the framework of fluid-mechanical systems (Toh 
& Itano 1999, Skufca et al. 2006, Schneider et al. 2007, Schneider et al. 2008, Duguet 
et al. 2007, Viswanath 2007, Viswanath & Cvitanovic 2008). The familiar concept of 
basin boundaries that separate different attractors was extended to the situation of 
a saddle coexisting with an attract or. We showed that the orbits defining the basin 
boundary are a set that changes smoothly when crossing a crisis point where one of 
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the attractors looses its stability. Beyond crisis we denote the closure of this set as the 
edge of chaos. The edge can be tracked by an iterative algorithm that exploits local 
properties only and hence can be used both in the situation of co-existing attractors as 
well as transient chaos coexisting with an attractor, see figure |6]( a, right panel). 

A standard procedure to determine the basin boundary is backward-iteration. It 
is more efficient than the direct forward sampling of phase space which was used to 
generate figure [S] The effort of backward iteration to determine a boundary of box- 
counting dimension Db with a resolution e scales like e~^^ . In contrast, the direct 
iteration scales quadratically with the resolution, i.e., like The edge tracking 

algorithm adopted in the present work (figure [g]) roughly requires the same numerical 
effort as backward iteration of the boundary of the region, and it has the additional 
benefit that beyond the crisis it focusses on the dynamically most relevant region of the 
edge of chaos while the backward iteration also tracks the circumference of all fingers 
shown in figure [7| 



6.2. Geometry of the boundary 

The geometry of the boundary separating laminar and turbulent dynamics can be 
studied in lifetime plots, where lifetime of initial conditions is either analysed for fixed 
parameters as a function of state-space coordinates, or by varying a parameter and a 
coordinate. 

For fixed parameters the separating boundary can be smooth or rough. The analysis 
in section 4.2 shows that roughness can be observed only if (a) the dynamics in the 
edge is chaotic, and (b) the Lyapunov exponent characterising the chaotic dynamics 
in the boundary is larger than the one in perpendicular direction. Roughness of the 
boundary hence is an indicator that there is a strong chaotic dynamics in the basin 
boundary. Since there is no a priori reason why the Lyapunov exponent pointing out 
of the separating boundary should be large, it will be interesting to identify a fiuid 
mechanical realization of rough basin boundaries. Ideally, the system should have a 
control parameter that infiuences the ratio of the Lyapunov exponents in the longitudinal 
and transverse directions. A good candidate might be Taylor- Couette fiow between 
independently rotating cylinders with a narrow gap, in which case it is close to the 
planar shear fiows mentioned earlier (Faisst & Eckhardt 2000). But it might also be 
possible to find evidence for rough boundaries in other parameter regions and geometries 
where a multitude of attractors can coexist (Abshagen et al. 2005). 

We have shown here how features of the boundary in the phase space relate to 
features in the parameter-coordinate space. The latter representation is typically studied 
in hydrodynamic systems where the Reynolds number Re is adopted as parameter. 
Increasing Re the boundary shows folded hierarchical organised tongue-like structures. 
In our model they appear shortly before or at the parameters of the boundary crisis 
of the turbulent attractor. The tongues have thus been related to the emergence of 
dynamical connections (c/. Rempel et al. 2004) between the relative attractor on the 
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edge of chaos and the attractor mimicking stable turbulent motion. These fingers result 
from the chaotic motion of the attractor undergoing a crisis. The presence of similar 
tongue-like structures in hnearly stable shear flows (Darbyshire & MuUin 1995, Moehhs 
et al. 2004a, Moehlis et al. 2004&, Schneider et al. 2007) further supports the idea of a 
turbulence generating chaotic saddle in these flows. The long persistence of turbulent 
motion, i.e.^ its tiny decay rate, may then be interpreted as another manifestation of 
supertransients (Lai & Winslow 1995, Breban & Nusse 2006). 

The local attractor embedded in the separating boundary - the edge state is an 
object both of theoretical and practical interest. The model shows that the local 
attractor can be a fixed point, a periodic orbit or a chaotic set. The type of dynamics 
in the boundary can be chosen independently of whether turbulent motion is generated 
by an attractor or a saddle. Thus, it is not a priori clear which type of edge state one 
should expect in transitional shear flows. A chaotic edge state has been identifled in pipe 
flow (Schneider et al. 2007), and a simple flxed point in plane Couette flow (Schneider 
et al. 2008). However, based on our present model we expect that other flow geometries 
show edge states with various other types of dynamics. 

6.3. Outlook 

The iterated edge tracking algorithm can be used to analyse any dynamical system 
showing two coexisting types of dynamics (Cassak et al. 2007). Without additional input 
the method can be used to analyse the position of the boundary and of trajectories in the 
boundary. A promising future application might be in control strategies, where the edge 
tracking is used to identify target states for chaos control (Schuster 1999). In various 
technological applications one is interested to intentionally induce turbulence or keep 
the flow laminar (Bewley et al. 2001, Hogberg et al. 2003, Kawahara 2005, Fransson 
et al. 2006, Wang et al. 2007). Up to now the setting up of the required effective 
control mechanisms mostly relies on empirical strategies, long-term experience and 
intuition. The edge tracking mechanism can provide additional guidance by identifying 
flow structures on which actuators could focus. 

6.4' Closing remarks 

The concept of the edge of chaos provides a powerful framework to analyse nonlinear 
dynamical systems where attractors coexist with a chaotic saddle and where the 
traditional concept of basin boundaries can no longer be applied. The approach still 
works for systems with several positive Lyapunov exponents. In that situation it 
provides insight into local attractors in the edge of chaos. 
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